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■ ABSTRACT 

This is the fourth paper in a series in which we attempt to put constraints on the flattening of dark halos in disk galaxies. We observed 
for this purpose the HI in edge-on galaxies, where it is in principle possible to measure the force field in the halo vertically and 
radially from gas layer flaring and rotation curve decomposition respectively. As reported in earlier papers in this series we have for 
this purpose analysed the HI channel maps to accurately measure all four functions that describe as a function of galactocentric radius 
the planar HI kinematics and 3D HI distribution of a galaxy: the radial HI surface density, the HI vertical thickness, the rotation curve 
and the HI velocity dispersion. In this paper we analyse these data for the edge-on galaxy UGC7321. 
We measured the stellar mass distribution (M = 3 x 10 8 M Q with M/L R £, 0.2), finding that the vertical force of the gas disk dominates 
f ) ' the stellar disk at all radii. Measurements of both the rotation curve and the vertical force field showed that the vertical force puts a 

much stronger constraint on the stellar mass-to-light ratio than rotation curve decomposition. Fitting of the vertical force field derived 
. from the flaring of the HI layer and HI velocity dispersion revealed that UGC7321 has a spherical halo density distribution with a 

flattening of q = c/a = 1.0 ± 0.1. However, the shape of the vertical force field showed that a non-singular isothermal halo was 
required, assuming a vertically isothermal HI velocity dispersion. A pseudo-isothermal halo and a gaseous disk with a declining HI 
velocity dispersion at high latitudes may also fit the vertical force field of UGC7321, but to date there is no observational evidence 
that the HI velocity dispersion declines away from the galactic plane. We compare the halo flattening of UGC7321 with other studies 
in the literature and discuss its implications. Our result is consistent with new n-body simulations which show that inclusion of 
hydrodynamical modelling produces more spherical halos. 
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ON 1 ■ Introduction galactocentric radius, which we used to measure the HI flaring 

■ of each galaxy in our sample. 
CO ; In paper I in this series dO'Brien et aljUoiOa) we presented HI i n the present paper we have fitted the vertical shape q = 
^ ! observations of a sample of 8 edge-on, HI rich, late-type galax- c / a f the halo density distribution for the northern galaxy 
O ' ies - The aim of the P ro J ect has been described there in detail. UGC7321. This particular system was chosen as a first ap- 
Q ; Briefly, we attempt to put constraints on the flattening of dark ha- plication since the sensitivity of the HI imaging obtained for 
. los around disk galaxies by measuring the force field of the halo UGC7321 at the VLA was 5 times greater than that for the 
•• ■ vertically from the flaring of the HI layer and radially from ro- southern galaxies that we observed with the ATCA, allowing 
. ^ ! tation curve decomposition. For the vertical force field we need more accurate measurement of the gas layer flaring to high lat- 
^ . to determine in these galaxies both the velocity dispersion of the i tu des, and better measurement of the HI velocity dispersion. 
5_h ; HI gas (preferably as a function of height from the central plane Due to its northern location, UGC7321 was not in our initial 
, f3 ; of the disk) and the thickness of the HI layer, all of this as a func- southern galaxy sample for which we measured near-IR and op- 
tion of galactocentric radius. In addition we also need to extract tical stellar photometry at Siding Spring O bservatory. However , 



information on the rotation of the galaxy and the deprojected HI VLA HI data observedby Lyn Matthews (Matth ews et all 1 1999 ) 
surface density, also as a function of galactocentric radius. was available for this galaxy and Michael Pohlen dPohlen etal 



In paper II (lO'Brien et al.L 1201 Obi) we discussed methods to 120021) kindly supplied fl-band photometry which allowed us to 

analyse the HI observations in edge-on galaxies and presented a derive the stellar luminosity density necessary to analyse the 

new method to measure the radial distributions, rotation curves halo density distribution. 

and velocity dispersions. We applied this metho d to our sam- In Sect.|2]we present the surface brightness and deprojected 

pie of galaxies in the third paper in this series: lO'Brien etalj luminosity volume density, and our derivation of the halo core 

( 1201 Pel) . In that paper we also developed a new method to derive radius, halo asymptotic velocity, and stellar mass-to-light ratio 

the thickness of the HI layer, or 'flaring profile', as a function of by rotation curve decomposition. In Sect. [3] we present a new 

simple method used to measure the halo shape using the ver- 

tical gradient of the vertical force, dK z /dz, with the usual as- 

* For correspondence contact Ken Freeman or Piet van der Kruit sumption of gas in hydrostatic equilibrium to determine the total 
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X (kpc) 

Fig. 1. Projected stellar surface density of UGC7321 averaged 
over all quadrants. Contours are separated on a log scale at 0.1, 
0.3, 1,3, 10, 3OL pc- 2 . 



dK z /dz of the galaxy. The resulting halo shape measurement for 
UGC7321 is presented and discussed in Sect. [4] and compared 
to other measurements of dark halo flattening in Sect. Sect. [6] 
summarizes our conclusions. 



2. Stellar surface brightness and deprojected 
luminosity density of UGC7321 

Fig.Q]shows the /?-band surface brightness - averaged over four 
quadrants - of UGC7321 with contours ranging from 0.1 to 30 
L pc -2 in steps of 0.5 dex. The observati ons and photometri c 
calibration o f UGC 7321 are discussed in Pohlenetal. (2002). 
iPohlen et al.l (l2003h analysed the projected surface brightness 
showing peanut-shaped deviations from elliptical fits to the 
isophotes at z heights greater than 0.5 kpc above the plane. These 
deviations provide strong evidence of a stellar bar, although it is 
difficult to measure the scale of the bar from the scale of the 
boxy-peanut shaped bulge (Athanassoula, private communica- 
tion). 

Using an exponential radial surface profile we fit the central 
surface brightness and apparent radial scale length of the pro- 
jected surface density. In logarithmic units 



fi(X) = fio + 1.086(X/7zx) mag arcsec 



-2 



(1) 



where X denotes the major axis distance. We found the pro- 
jected central /?-band surface brightness to be 22.0 mag arcsec -2 , 
with a scale length of 4.0 + 0.3 kpc, in agreement with that 
measured by Pohlen et al. (2003). The projected central surface 
brightness is consisten t with the B and 7?-band measurement by 
Matthews et al] ( 1 19991) when the internal extinction model that 
Matt hews et alJ used is factored in. The projected surface bright- 



ness shows a small nuclear feature smaller than 1 kpc, and an 
exponential profile from 0.5 to 6.5 kpc, declining steeply after 7 
kpc. 

2.1. Deprojection 

To deproject the luminosity distribution from the edge-on pro- 
jection, we assume azimuthal symmetry and perform a direct 
deprojection of the projected surface density on the sky using 
the inverse Abel transform 



I(R,z) = — 

n 



1 



dI(X,z) 



V* 2 - R 2 dX 



dX, 



(2) 



where X is the position along the major axis and R is the galac- 
tocentric radius in the cylindrical coordinate system. Applying 
the inverse Abel tranform to each z-plane yields the luminous 
volume density I(R,z) of the galaxy in L pc -3 as a function 




Fig. 2. Stellar luminosity volume density in R and z after de- 
projecting the observed surface density using the Inverse Abel 
Integral. Contours are separated on a log scale at 0.0001, 0.003, 
0.001, 0.003, 0.01, 0.03, 0.1, 0.3 L pc -3 . The radial and vertical 
scalelengths, as shown in Fig. [3] are Iir = 2.65 kpc and h z - 245 
pc. 



of cylindrical radius R and height z. To minimise the incidence 
of local regions with a negative derivative dI(X,z)/dX, we per- 
formed the deprojection after smoothing the projected surface 
density by several different 2D Gaussians with FWHM ranging 
from 256 pc to 1027 pc (or 10 and 40 pixels, respectively). 

Fig.|2]shows the derived volume density I{R, z) of UGC7321, 
while Fig.|3]shows the radial and vertical surface brightness ob- 
tained by integrating over the volume density. UGC7321 is in- 
deed a very low surface brightness galaxy. The face-on radial 
scalelength is Iir = 2.65 ± 0.17 kpc. Both the volume density 
and the radial profile, show that UGC7321 clearly has a small 
central nuclear region that is approximately 8 times brighter in 
surface brightness than the fitted central surface brightness of 
Hr = 23.4 + 0.14 mag arcsec -2 derived assuming an exponen- 
tial disk only. The peak face-on central surface brightness of the 
nuclear region is /jr =21.1 mag arcsec -2 . 

The central luminosity volume density of UGC7321 is 0.3 
L pc -3 in the central nucleus averaged over 250 pc. At 2.2 scale- 
lengths the midplane volume density is 0.002 L pc -3 , much less 
than the luminosity density near the Sun 0. 1 L pc -2 ) as would 
be expected for such a low surface brightness galaxy. 



2.2. Rotation curve decomposition 

In the standard manner we decomposed the rotation curve to ob- 
tain the parameters of a spherical pseudo-isothermal halo, and 
constrain the stellar mass-to-light ratio. The radial surface den- 
sity of the stars and gas was used directly to derive the rotation 
curve contribution due to each luminous mass component, and 
the observed rotation curve fitted such that 



v,, = v; 



»/>. 



(3) 



where v b s (R) is the observed rotation curve, and vi,(R), v s (R) 
and v g {R) are the rotation due to the halo, stellar and gas mass 
components, respectively. The radial surface density of each lu- 
minous component was obtained by integrating over z, and a 
constant /?-band stellar M/Lr was adopted. The HI surface den- 
sity was scaled by a factor of 1 .4 to include He and (a minimal 
amount of) molecular gas. 

The observed rotation curve was fitted using a spherical 
pseudo-isothermal halo density distribution 



Ph(R) = Ph,o 



R 



R 2 +R 2 



(4) 
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Fig. 3. Exponential fits to the radial and vertical surface bright- 
ness profiles formed by integrating over the deprojected volume 
density. The arrows mark the radial scalelength (left) and vertical 
scaleheight (right). 



with corresponding rotation curve 



v 2 h (R) = AnGp K0 Rl 



i R < t I R 
1 arctan — 

R \R C 



(5) 



The pseudo-isothermal halo shown above is defined by the core 
radius R c and central density and asymptotes at R — > oo to 



v h,oo 



= AjiGp hfi R 2 c , 



such that the rotation can also be written 



1 R * * I s 
1 arctan — 



(6) 



(7) 



By definition the rotation curve measures the total force in 
the radial direction Kr = v 2 (R)/R; however it is unable to con- 
strain the halo flattening q. 

By evaluating the radial force Kr(R, z) of a flattened pseudo- 
isothermal halo with density distribution 



P(R,z) 



Ph,oR c 



(8) 



R 2 C +R 2 +Z 2 /q 2 ' 

ISackett et alj (Q994) show that the corresponding rotation curve 
in the midplane is 



r 



arctan y \R 2 + a 



q R c 



arctan 



R 2 



q 2 R 2 



1/2 



where 



y = 



a = (l-q 2 )R 2 



(9) 



(10) 



and the asymptotic rotation M (^) is very similar to the asymp- 
totic rotation of a spherical isothermal halo with 



h,oo 



(q) = 4nGp h , R 2 c f(q), 



where 
/(<?) = 



arccos(g). 



(11) 



(12) 



We use the rotation curve of a spherical pseudo-isothermal 
halo, as the shape of the rotation of a similar, but flattened halo 




RN/r 



Fig. 4. Scaled rotation curves of a flattening pseudo-isothermal 
halo for two different values of the halo flattening q. The thick, 
black curve shows the rotation curve of q — 0.3 and the thin, 
red curve shows the rotation for q = 0.9 with the radius scaled 
by 0.84. The vertical axis shows the rotation normalised by the 
asymptotic rotation vi h00 (q), while the abscissa shows the radius 
normalised by the core radius R c . The nearly identical shape of 
the two curves shows that shape of the rotation curve of the halo 
is almost independent of q, with the radial scaling varying by a 
constant factor (of only * 15% in this case) over a large range of 
q. 



is almost the same. This similarity implies that the measured 
asymptotic rotation derived from a spherical pseudo-isothermal 
fit to the rotation curve also defines the asymptotic rotation of a 
flattened pseudo-isothermal halo via Eqn. (fTTl i. 

In Fig. [4] we show how the rotation curve of a flattened 
pseudo-isothermal halo varies with q (for q < 1) at z=0. The 
vertical axis shows the rotation normalised by the asymptotic ro- 
tation v/, oo(g), while the abscissa shows the radius normalised by 
the core radius R c . The thick, black curve shows the rotation for 
q = 0.3, while the thin, red curve shows the rotation for q = 0.9 
with the radius scaled by 0.84. The nearly identical shape of the 
two curves shows that shape of the rotation curve of the halo is 
almost independent of q, with the radial scaling varying by only 
=s 15% over a large range of q. The derived halo central den- 
sity does vary significantly with q, becoming denser for more 
flattened halos: 



Jq - 1) i 

AnG R 2 f( q y 



(13) 



In Sect. |4] we measure q by fitting the gradient of the ver- 
tical force dK-Jdz, where the gradient of the flatte ned pseudo- 
isothermal halo is also given bv lSackett et all ( 1 1994b . As K Z (R, z) 
also depends on the same parameters q, halo core radius R c and 
halo central density p/, o, we can use the asymptotic halo rotation 
Vh,oo and the core radius R c which are well measured from rota- 
tion curve decomposition, and just perform a single parameter 
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fit to the flattening q, as the central density p/, o is tied to Vh,<x>, R c 
and q. 

In Fig. [5] we show the best fitting stellar and halo rotation 
curves derived from rotation curve decomposition of UGC7321. 
The observed rotation curve is shown by the thick black line, 
while the rotation due to gas, stars and the halo are shown by 
dashed lines (from the bottom up stars (red), gas (green) and 
halo (blue). The resulting fit is shown by the grey full-drawn 
line (yellow). The best fit was achieved with a sub-maximal (see 
below) stellar M/Lr of 1.05, which scaled the stellar component 
from a luminosity of L = 4.0 x 10 8 L to a mass of 4.2 x 10 8 
M . UGC7321 is a gas-rich low surface brightness galaxy with 
a Mhi/Lr = 2.2; thus even with an M/Lr near unity, the stellar 
mass is roughly a third of the gas mass (M g = I AM hi)- 

The plot shows that the rotation curve is well fit by a pseudo- 
isothermal halo with core radius R c = 0.52 ± 0.02 kpc, and cen- 
tral density p/,,o = 0.73 ± 0.05 M pc~ 3 . As the halo core radius 
and halo asymptotic rotation are nearly independent of the halo 
shape in flattened pseudo-isothermal halos, we adopt those mea- 
surements from rotation curve decomposition for our halo shape 
modelling in Sect. [3] The halo is clearly the dominant component 
at all radii, comprising over 90% of th e total galaxy mass o f 3.2 X 
10 10 M Q at the last measured point. lUson & Matthewi d2003l) 
reach a very similar conclusion. UGC7321 has a similar dark-to- 
light mass ratio, although it is significantly more m assive than 
DDQ 154 and shows evidence of a bar dCarignan & Freemanl 
1988). 

The dominance of the halo at all radii means that the stellar 
mass-to-light ratio is rather unconstrained. Forcing a maximum 
disk fit yielded a stellar M/Lr of 2.9, but resulted in a very poor 
fit with rotation velocities in error by 15-25% in the inner disk. 
Even with a maximum disk fit, the peak rotation for the stellar 
disk was only 43 km s _1 . The best fit occurs for an M/Lr of 1.05, 
but all M/Lr values less than 1.05 provide acceptable fits to the 
rotation curve. Consequently, we use 1 .05 as an upper limit in 
the analysis in Sect. [3] 

The low mass-to-light ratio and extremely sub-maximal na- 
ture of the UGC7321 stellar disk implies that the recent star 
formation rate excee ds the average rate ove r the galaxy's life- 
time. Observations by Matthews et al. (1999) that found a signif- 
icant fraction of young stars support our results. Detailed studies 
of the vertical dis k structure indicate multiple disk components 
(Matthews, 2000). Our finding of a very low mass-to-light ra- 
tio (in the /?-band) warrants further studies of the structure and 
composition of the stellar disk, but this is beyond the scope of 
the present paper. 

3. Method for fitting the halo shape 

The thickness of the gas layer depends on the vertical force K 7 , 
and hence on the shape of the d ark matter hal o (iMalonevi 1 1 992 ; 
Oiling & van Gorkoml 1 1993b iKundic et all Il993t IMalonevi 
19931) . A gas layer with less flaring for a given gas velocity dis- 
persion implies a stronger K z . 

Analytically it has been shown that at large radii the thick- 
ness of the gas layer is roughly proport i onal to the sq uare root 
of the halo flattening q (Malonev, f 19931: lOllingl 119951) and that 
the flaring should i ncrease radially in an exponential fashion 
dvanderKruitlll98lh . This was confirmed in our measurements 
and the earlier study of lRupenl d!991l) comprising high resolution 
VLA HI observations of NGC891 and NGC4565. 

We determine the halo flattening by measuring the z-gradient 
of the total vertical force K z tot (R, z) from the equation of hydro- 
static equilibrium for the gas layer, and evaluating dK z /dz for 
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Fig. 5. Rotation curve decomposition of UGC7321. The ob- 
served rotation curve is shown by the thick (black) line, while 
the rotation due to the other components are shown by dashed 
lines; from the bottom up stars (red), gas (green) and the halo 
(blue). The disk has in this fit an M/Lr of 1.05. The resulting fit, 
vK it = Vg + v 2 s + vjj, is shown by the grey full-drawn line (yellow). 



each luminous mass component using Poisson's equation. The 
halo gradient dK z j,/dz is modelled using the equation for the 
vertical force K 7 of a flattened pseudo-isothermal halo given in 
ISackettetal.ldl994l) . 

Given the gas disk is in equilibrium, the gas pressure gradient 
and internal forces must exactly balance the gradient of the total 
gravitational potential of the galaxy, where the total gravitational 
potential <t>, , is the sum of the stellar, gas and halo potentials, 
®tot = + ®s + ®g- Assuming that the gas velocity dispersion 
is isothermal in z (though not in R), the equilibrium condition in 
z becomes 



-^[al g (R)lnp g (R,z)] = 



dK z (R,z) 
dz ■ 



(14) 



If we further assume that the gas density distribution is 
Gaussian in z, the vertical gradient of the total K z becomes a 
simple function of the gas velocity dispersion and the vertical 
FWHM thickness of the gas, both functions of radius which we 
measured in paper III, namely 



dK ZJO ,(R,z) 
dz 



K g (R) 



(15) 



(FWHM z>g (tf)/2.355) 2 ' 

From this we see that the vertical gradient of K ZJOt derived from 
such a gas disk is constant in z, and varying in R. 

The gradient of the vertical force of each of the luminous 
components was directly calculated from the Poisson equation 
for each component 

dK z (R,z) ,_^_,„ _ x . 1 d(RK R ) 



dz 



-4nGp(R,z)- 



R dR 



(16) 



where p(R,z), Kr, and K z correspond to the volume density and 
forces of that component. In terms of v 2 (R), the squared rotation 
due to that mass component, we can rewrite this as 

-AnG P {R,z) + i^. (17) 



dK z (R,z) 



dz 



R dR 
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Fig. 6. The vertical gradient of the vertical force. The total gradient, dK zwt jdz determined from hydrostatics is shown as the thick, 
black line (labelled 'total'), while the gradients for the gas and stars is shown in green and red (and labelled 'gas' and 'stars'), 
respectively. The stellar gradient dK z ^/dz is calculated using the stellar mass-to-light ratio of M/Lr = 1.05 ± 0.33, derived from 
rotation curve decomposition. This value of the stellar M/Lr is clearly too high as the gradient allowed for the halo (shown in 
blue and labelled 'halo') is positive at all radii, requiring negative halo densities. The stellar M/Lr needs to be reduced to allow a 
negative halo gradient dK z j,jdz. The cyan line (labelled 'halo+stars') shows the sum of stellar and halo components determined by 
dK z j t/dz - dK zg /dz, which must be fitted with a non-negligible stellar M/Lr and a positive density halo. 



Consequently, the gradient of the halo force must satisfy 



dK zJl (R,z) dK Uot (R,z) dK z , s (R,z) 



dz 



dz 

dK z JR,z) 
dz ' 



dz 



(18) 



where the subscripts s, g and h, denote the stars, gas and halo, 
respectively. 

The vertical dK z j,(R, z )/dz force of the fl attened pseudo- 
isothermal halo is given in lSackett et alj d 19941) . As the asymp- 
totic halo rotation and the halo core radius were well determined 
from the rotation curve decomposition, and the central density 
determined by q, the fitting of dK Z: i,/dz reduces to a fit with a 
single parameter q. 

Comparison of the stellar vertical force gradient with the to- 
tal vertical force gradient shows that the vertical force puts a 
much stronger constraint on the stellar mass-to-light ratio than 
does the radial force fitting undertaken in rotation curve decom- 
position. Inclusion of the gas self-gravity requires that the stellar 
vertical force gradient dK z s (R, z)/dz must be 



dKUz) dK Uot {z) dK Ztg (z) 



dz 



dz 



dz 



(19) 



where from here on we take the ^-dependence of K z as implicit 
and write the derivatives as total derivatives. 

Given that the stellar mass density and its squared rotation 
are both proportional to the stellar mass-to-light ratio, we see 
from Eqn. ( ITTb that the vertical gradient of the vertical force of 



the stars is linearly related to the stellar mass-to-light ratio. As 
the low stellar luminosity meant the rotation curve decomposi- 
tion was relatively insensitive to the stellar mass-to-light ratio, 
we consider the stellar mass-to-light ratio to be a free parameter 
and fitted 



dK 



Z.IOI 



dK 



dz dz 

as determined from the observations as a function of R with 



dK 



dz 



-./, + d K : s 



dz 



where the K Zt h is modelled by the flattened pseudo-isothermal 
halo and the z-gradient of the stellar vertical force is 



dK z , s (M\ dK zl 



dz 



z \L 



dz 



M/L,< = 1 



The total K z gradient and the gas K z gradient come directly from 
observations without any free parameters. 

We recall from above that the vertical gradient of K zJot de- 
rived for a Gaussian gas disk is independent of z. To be most sen- 
sitive to the constraints from the luminous mass density, we fit 
the gradient of the vertical force near the midplane at z = 100 pc, 
high enough to avoid the bulk of the internal extinction caused 
by dust in the plane of the thin stellar disk. 
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4. The halo shape of UGC7321 

4.1. Results of the fitting 

Fig. [6] shows the vertical gradient of each of the K z components 
for UGC7321. The thick, black line (labelled 'total') is the gra- 
dient determined from the gas flaring and velocity distribution 
using Eqn. HSi . The gradient due to the stellar disk is shown by 
the red line (labelled 'stars') and has been calculated with the 
mass-to-light ratio of 1.05, which was the best fitting value in 
our rotation curve decomposition. We see immediately that the 
stellar dK zs jdz alone is comparable to the total dK ZJOt /dz given 
by the hydrostatics over most of the range of R. The gradient due 
to the HI is the line labelled 'gas' (green). We can subtract this 
gradient for the gas from the total gradient and derive the gradi- 
ent due to the sum of the halo and the stars. This is the (cyan) line 
labelled 'halo+stars'. Subtracting the gradient for the stars from 
the gradient for the (halo + stars) then leaves the gradient that 
should be attributed to the halo alone (blue line labelled 'halo'). 
With the adopted M/L ratio for the stellar disk, the halo gradient 
turns out to be positive, which is unphysical. 

From this example, it is clear that the gradient of K ZJOt mea- 
sured from the hydrostatics provides a very strong constraint on 
the stellar M/Lr ratio. Even with a zero-mass halo, which is ex- 
cluded by the rotation curve fit, we see that the stellar M/Lr must 
be less than 1 to leave room for the gradient of K z given by the 
gas self-gravity. With the necessary inclusion of the gas, we find 
that the stellar mass-to-light ratio M/Lr must be i 0.6. The dif- 
ference dK ZJO tldz - dKg/dz (cyan; labelled 'halo+stars') consti- 
tutes the combined dK z /dz of the halo (blue; labelled 'halo') and 
stars (red; labelled 'stars'). The rotation curve decomposition re- 
quires a positive halo mass density at all radii, thus dK Zt h/dz < 
for all radii. This constraint requires that the stellar M/Lr <k 1 . 

Inspection of the dK zto ,/dz - dK Ztg /dz difference (cyan; 
'halo+stars') shows a steep gradient at small R, flattening at large 
radii, particularly where the gas layer undergoes exponential 
flaring at radii outside 7 kpc. We see from Fig.|7]that this is sim- 
ilar to the characteristic shape of dK z j,/dz for pseudo-isothermal 
halos with different halo flattenings. But despite this similarity, it 
was not possible to fit the halo flattening q while holding the core 
radius R c and the asymptotic halo rotation v/ !jCO fixed, even with a 
zero mass stellar distribution (stellar M/Lr = 0) and allowing q 
to range between oblate and prolate shapes. By adjusting q, and 
keeping M/L#| storj small, it is possible to get a similar shape to 
the difference gradient (cyan; 'halo+stars'), but it is always off- 
set to larger negative values of the gradient. This implies that the 
asymptotic halo rotation scale derived from the rotation curve is 
too high, as the magnitude of dK z j,/dz is proportional to the v 2 h M 
(see ISackett etallll 994 . 

This may be an artefact of the adopted pseudo-isothermal 
halo model. While the flattened pseudo-isothermal (PIT) model 
is computationally convenient, we note that a true (spherical) 
non-singular isothermal (IT) model was initially adopted for 
dark halos (Carignan & Freeman 1985, 1988). Few studies have 
compared the relative merits of the PIT and IT dark halo mod- 
els. In their paper on halo s caling laws for disk galaxies (S c and 
later) and dwarf speroidals. iKormendv & Freemanl (12004) com- 
pare halo fits to rotation curves over a large sample and generate 
scaling laws between halo parameters measured with a IT halo 
and those with a PIT halo. In Fig. [8] (adapted from Fig. 1 of 
Kormendy & Freeman, 2004), we show IT and PIT halo rotation 
scaled to the same asymptotic rotation. As can be seen, the rota- 
tion of the IT halo rises above the asymptotic rotation speed be- 
fore declining to it at large radii, while the rotation curve of a PIT 
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Fig. 7. The vertical gradient equated near the midplane of a flat- 
tened pseudo-isothe rmal halo, as calculated from Eqn. (6) of 
Sackett et al. (1994) for both oblate and prolate halos ranging 
from q of 0.5 to 1 .5 in steps of 0. 1 . The bottom curve corresponds 
to the most flattened halo with q - 0.5, becoming increasingly 
shallow as the halo gets less flat. The gradient of a spherical 
halo is shown by the dashed line. The top curve corresponds to 
2 = 1.5. 



halo approaches the asymptotic rotation from below. The declin- 
ing shape of the rotation curve for the IT model would provide 
a lower and possibly more realistic estimate of the asymptotic 
rotational velocity v/, jC „ from a rotation curve decomposition of 
rotation data which in practice does not extend to radii in excess 
of a few halo core radii. 

If an IT halo was fitted to the observed rotation curve of 
UGC7321, the asymptotic rotation would be approximately 20- 
40% lower than that of the PIT halo. This would provide the 
lower asymptotic rotation scale v/ 100 necessary to fit the dif- 
ference dK-Jdz of UGC7321 with a flattened halo over radial 
ranges from 1.5 to 9 kpc. 

Flattened non-singular isothermal halos could be formed by 
halo rotation or anisotropy of the velocity dispersion. The rota- 
tion is unlikely to be figure rotation, as figure rotation of triaxial 
halos measured in n -body simulations was found to be very slow 
dBailin & Steinmetzl 12004]) (0.148/j km s" 1 kpc 1 , where h is 
//o/lOO, insufficient to flatten halos more than q ~ 0.7). The ve- 
locity dispersion anisotropy of the halo dark matter would allow 
either prolate or oblate halos, just as velocity anisotropy of the 
stars in the brighter elliptical galaxies defines the galaxy shape. 

With the asymptotic rotation as a free parameter in addition 
to the stellar M/Lr and q, we found that the residual dK z /dz 
curve ('halo+stars', cyan) is best modelled with a halo shape of 
q = 1.0 + 0.1. Robust least squares minimization fitting using a 
Levenberg-Marquardt algorithm (MINPACK- 1 ) favoured a zero 
mass stellar disk, but fits were almost as good for an M/Lr = 0.2 
stellar disk. These fits were successful over the radial range from 
2-9 kpc. 
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Fig. 8. This figure shows the difference in the rotation curve and 
velocity dispersion of isothermal and pseudo-isothermal halo 
models. The top panel shows that rotation of an IT halo rises 
above the asymptotic rotation speed at large radii, before declin- 
ing to it, while the rotation curve of a PIT halo approaches the 
asymptotic rotation from below. Both rotation curves are scaled 
to the same asymptotic rotation speed to allow comparison. The 
lower panel shows that an isotropic pseudo-isothermal halo is 
not isothermal, exhibiting a lower velocity dispersion at low 
radii. (Adapted from Fig. 1 of Kormendy & Freeman, 2004.) 



We illustrate this first for the unphysical case where there is 
no mass in stars in Fig. [9] This figure is organised in the same 
manner as Fig. [6] The gradient due to the stars is now zero at all 
Z. Recall that the (cyan) line 'halo+stars' is the observed gradi- 
ent, which has to be fit. The smooth (also cyan) line 'halo' is that 
fit (the dashed -blue- line superimposed is that of the halo alone, 
which is the same when M/Lr is zero). This best fit was achieved 
with an asymptotic PIT halo rotation reduced by 30 ± 5% com- 
pared to the PIT fit to the rotation curve in Fig. [5] 

It is remarkable that the shape of the /^-dependence of 
dK z /dz for the adopted halo model in Fig.[9]agrees so well with 
the shape of the K- gradient derived from the HI flaring and 
velocity dispersion, at least for radii < 9 kpc. Although some 
rescaling of the strength of the K z force was needed, we see that 
the density distribution of the adopted spherical PIT, using the 
core radius derived from the rotation curve fit, also provides the 
correct radial variation of the K z gradient. This need not have 
happened. Although the Kr estimate from the rotation curve 
and the hydrostatic estimate of dK z /dz come from analysing the 
same XV data (see paper III), the two functions come from dif- 
ferent features in the XV data, so are relatively independent. 

In Fig. [10] we show the fit for a stellar M/Lr = 0.2 disk. 
The lines labelled 'halo+stars' (cyan) show the gradient as de- 
duced from the observed total gradient minus that of the gas in 
the full-drawn line (which of course is the same as in Fig. [9} 
and that of the sum of the gradient of the halo model fit and 
that deduced from the stellar distribution with M/Lr = 0.2 as 



the dashed line. The gradient from the stellar disk alone is the 
(red) curve labelled 'stars' and for the halo the (blue) curve la- 
belled 'halo'. For this case an asymptotic halo rotation reduc- 
tion of 50 ± 5% was needed. In effect, reducing v/, j0O and reduc- 
ing the stellar M/Lr have similar effects of increasing the mag- 
nitude of the asymptotic value of the difference dKJdz curve 
('halo+stars', cyan). For both cases (M/Lr — in Fig. [9] and 
M/Lr = 0.2 in Fig. |T0j, the shape of this difference dK z /dz 
curve dictated (given the derived core radius for the dark halo 
from our rotaton curve decomposition) a halo flattening close to 
spherical. 

At radii larger than 9 kpc, the strong flaring causes the dif- 
ference dK z /dz (cyan, 'halo+stars') to be too small to be fit with 
the same asymptotic halo rotation. Even the shallow gradient 
of dK z /dz given by a highly prolate halo, combined with a low 
Vh, w did not produce a good fit at these radii. We briefly discuss 
why the derived gradient for halo+stars may have been underes- 
timated. 

We have argued that a true isothermal halo may provide a 
more valid model for this analysis. Another possibility is that 
the gas velocity dispersion is not vertically isothermal. We were 
forced to adopted this assumption in the hydrostatic equation, 
because there is currently no available measurements of the z- 
dependence of the gas velocity dispersionQ Prior to our work the 
gas velocity dispersion had only been measured in a few face-on 
galaxies. Our HI disk modelling of edge-on galaxies has more 
than doubled the number of galaxies with radial gas velocity dis- 
persion measurements. These high resolution observations show 
that the gas velocity dispersion is not isothermal in radius, but 
its vertical properties are unknown and we had to assume that is 
is vertically isothermal. 

The gas velocity dispersion in disk galaxies is often ascribed 
to local heating by s upernovae and stellar winds in s tar forma- 
tion regions. Indeed. IShostak & van der Kr uit (1984) found in 
NGC628 that the gas velocity dispersion is systematically higher 
in the spiral arms than in between. On the other hand, similar ve- 
locity dispersions are seen in regions of star forma tion and in re- 
gions where there is no visible star formation (e.g. iMeurer et all 
1996), and both low and high surface brightness galaxies seem to 
have similar gas velocity dispersion. Sellwood & Balbus ( 1999) 
offer a plausible alternative, suggesting that weak magnetic 
fields in galaxies allow energy to be extracted from differential 
rotation via MHD-driven turbulence. This would result in a gas 
velocity dispersion that was proportional to the rotational shear 
due to the disk, resulting in similar gas velocity dispersion for 
galaxies with similar rotation curves. However, while heating 
caused by gas shear could generate radial variation in the gas 
velocity dispersion, it is unclear how it could cause the gas ve- 
locity dispersion not to be vertically isothermal. Conversely, the 
decline in star formation away from the midplane could cause a 
fall-off in gas velocity dispersion with z. 

A non-isothermal vertical gas velocity dispersion would 
probably have more of an effect at larger radii where the flar- 
ing gas is probing a larger range in z. A gas velocity dispersion 
declining with z would increase the absolute total vertical gradi- 
ent of K z derived from the equation of hydrostatic equilibrium. 
A significant increase of dK z lo Jdz would enable a larger asymp- 
totic halo rotation more consistent with a pseudo-isothermal halo 
and a larger stellar M/Lr. 



1 Because of S/N limitations, our measurement of the gas velocity 
dispersion in paper III models the Hi XV diagram integrated over z. 
Thus it represents a luminosity-weighted average dispersion as a func- 
tion of radius. 
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Fig. 9. The z-gradient of the vertical force showing successful fits of the halo shape for a zero mass disk M/Lr = 0. The black curve, 
labelled 'total' shows the total gradient dK zlot jdz derived from hydrostatics. The green ('gas') and red (coincident with the abcissa 
at the top of the figure) curves show the z-gradient derived using Poisson's equation for the gas and stellar mass components. The 
cyan ('halo+stars') line shows the difference dK ZJOt j dz - dK~ ug j dz- This line is fitted by the halo gradient dK z j,/dz which is modelled 
by a flattened pseudo-isothermal halo. To achieve a successful fit to the cyan line, the asymptotic halo rotation was scaled down by 
30%. The measured halo flattening was q = 1.0 + 0.1. The vertical dotted (cyan) lines show the radial regime used for the fit. 



Another, less plausible, explanation of the dK z /dz fitting 
problem at large radius, is that the gas-to-Hl ratio used to scale 
the HI density to account for He and H2 is not constant. This is 
unlikely as the He content is mainly primordial and well known 
from big bang nucleo-synthesis. As He accounts for 0.34 of 
the additional 0.4 fraction, it is unlikely that a radially declin- 
ing molecular hydrogen distribution could significantly reduce 
dK Zyg ldz thus allowing a higher difference dK z /dz (cyan). 



5. Comparison to other work 

We first review earlier work on the flattening of dark halos in 
spiral galaxies. The earliest concern was whether the dark matter 
indicated by fiat rotation curves resided indeed in a more or less 
round halo or was part of the disk. That the latter was not the 
case was shown in 1981, using evidence from bu lge isophotes 
in ext ernal galaxies and star counts in our Galaxy ( Monet et al.L 
[1981 and from HI flaring in NGC 891 dvan der KruiA Il981h . 
Next, the question of the actual flattening q — cja of dark halos 
in spiral galaxies arose and we will now review previous work on 
this subject, starting with our Galaxy. One of the early methods 
is the analysis of the local surface density in the Solar neig hbour- 
hood using stellar kinemati cs. With this metho d iBinnet et alj 
d 19871) find 0.3 ^ q < 0.6, Ivan der MareH d!991l) q > 0.34 and 
iBienavme et"ai1 d2006h q > 0.5. 

At large radial distances of 8.0 < R < 60 kpc, RR Lyrae 
stars show the dark matter distribution to be flattened by q ~ 
0.7 dAmendt & Cuddefordlll99 4). Hyper-velocity stars open an- 
other promising way of probing the shape of the Galactic dark 
matter distribution. One star, assuming it is 70 kpc away, gives 



0.5 < q < 1.6 dGnedin et all 120051) . ISamurovic et ail if 1999) 
used the microlensing optical depth towards the Galactic bulge, 
LMC, SMC and M31 to probe the shape of the Galactic halo to 
large radii (R < 5/?2s)- However, they were not able to derive 
strong constraints: q = 0.6 ± 0.4. 

Since the discovery of the Sagittarius dwarf galaxy, mod- 
elling of its extended stellar tidal deb ris stream has b ecome 
one of the most promising methods. iMaiewski et al] d2003h 
show that the Sagittarius stream traces a great circle around our 
Galaxy, extending to radii of 2/?25 from the Galactic centre. If 
the tidal debris has made several orbits, the Galactic halo must 
be near-spheric al so that the strea m does not precess away from 
a single plane. iMerrifieldl d2004) argued that the apparent co- 
herency of the carbon star kinematics in the stream suggest that 
all the stars are on the same wrap, m aking it impossib le to con- 
strain the halo flattening. Converselv. llbata et al.l d2001l) contend 
that the stream has made several orbits, and from this infer that 
the Galactic halo must have flattening q > 0.7 in the radial range 
16 < R < 60 kpc. 

Recently, nu merical l y mod elling of small satellite infall on a 
Sgr-like orbit bv lHelmil d2004al) finds that tidal streams younger 
than about 2 Gyr lead to spatially coherent streams for a large 
range o f halo flattenings 0.6 < q < 1.6. Since then she (Helmi, 
l2004bh has significantly revised her initial measurement to a 
highly prolate shape with 1.25 i q i 1.5 b y constraining th e 
star samp le to the older S gr str eam stars of lLaw et al] (|2005). 
However, iJohnston et al.l ( 20051) dispute this result, finding a 
near- spherical halo with q ~ 0.83 -0.92. In a more recent analy- 
sis of the Magellanic Stream Riizicka et al. (2007) find a flatten- 
ing of 0.74 < q < 1.20. 
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Fig. 10. The z-gradient of the vertical force showing successful fits of the halo shape for stellar mass-to-light ratio M/Lr = 0.2. 
The curve are coloured as in Fig. [9] except that the line for the stars is now at non-zero values. The cyan line, labelled 'halo+stars', 
shows again the difference dK ztot /dz - dK zg /dz which is fitted by the stellar gradient dK zs /dz derived from Poisson's equation 
and the halo gradient dK z j,/dz which is modelled by a flattened pseudo-isothermal halo. To achieve a successful fit to this line, the 
asymptotic halo rotation was scaled down by 50%. The thin full-drawn (blue) line, labelled 'halo' corresponds to this halo, while 
the dashed (cyan) line shows the gradient due to this halo and that of the stars together. A stellar mass-to-light ratio significantly 
larger than 0.2 was not possible. Again, the measured halo flattening was q = 1.0 ± 0.1. The vertical dotted (cyan) lines show the 
radial regime used for each fit. 



The situation for halo shape measurement in external galax- 
ies is just as confusing, because some methods are suited only 
to specific types of galaxies. The determination of halo shape 
from polar ring galaxies is such a case. By simply comparing the 
equatorial and polar rotation curves it is possible to ascertain the 
flattening of the total potential. Using this method NGC4650A 
and AO 136-08 01 were found to be moderately flatten ed with 
q I 0.7 dWhitmore et all 119871; ISchweizer etail 119831 respec- 
tivel y), while MCG-5-7-1 w as found to be approximately spher- 
ical (Whitmore et al 1 119871) . A potentially more accurate method 
is to model the rotation along both axes using a multi com- 
ponent mass model comprising bulge, equatorial stellar and 
gas disks, and polar ste llar and gas rings. Using this method, 
ISackett & Sparkd (1 1990b originally found the halo flattening of 
NGC4650A to be 0.3 < q < 0.7; subsequent higher qual- 
ity observations w ere able to constrain the halo more tightly, 
to 0.3 < q < 0.4 ( iSackett et all 11994 . This method has also 
been applied t o AM2020-504 , wher e the flattening was found 
to be q ~ 0.6 CArnaboldi et all fl993h . Another method involves 
modelling of the twisting c aused by precession of the ring. With 
some specific assumptions, Steiman-Camer on et al.l (11992) con- 
strain the flattening of the NGC4753 halo to be 0.84 < q < 0.99. 
Finally, using the twisting of the morphological minor axis of 
the disk plane away from the kinematic minor axis to model the 
velocity field of polar rings, the flatte ning of the dark halo of 
A0136-0801 was found to be q ~ 0.6 (ISackett & Poggdll995l) . 

Another method that has been used to measure halo flat- 
tening is strong gravitational lensing. An early study of a 



doub le lens sys te m co mprising two spirals found q > 0.4 
dKoopmans et all 119981) . More recently there have been 
two studies of multiple quad lens syste ms finding q ^ 0.4 
dRusin & TegmanL l2110h and q ~ 0.7 (ICohn & Kochanekt 
2004 ), and another an alysis of a double lens system 0.6 < q < 
0.7 dChae etaill2002l) . 

Warps in stellar disks (e.g. Reshetni kov & Combesl 1 19981) 
offer several mechanisms to probe the halo shape of spiral galax- 
ies. One method uses the precession of the warped disk to con- 
strain the halo flattening. It has been ap plied to NGC2903, yield - 
ing a halo flattening of q = 0.80±0.15 (Hofner & Sparkelll994l) . 

It is also possible to measure the mean shape of vast num- 
bers of galaxies via weak gravitational lensing. Measurements of 
about 10 s lensed syste ms against about 10 6 background galaxies 
dHoekstra et all 120051) find a mean projected halo ellipticity of 
0.20^ and a mean projected halo flattening of (q) = 0.66+^ 
(1-cr error) . However, a larger investigation of about 2 million 
lensed galaxies against 32 million background galaxies from 
the SDSS datas et found no strong eviden ce of flattening, with 
(q) = 0.99+^ (iMandelbaum et alll2006l) . 

The results of halo flattening studies so far do not reveal a 
consistent picture. We believe that the method of the flaring of 
the gas layer is among the most promising, at least f or late- 
type spiral galaxies. First tried by ICelnik et al.l d 19791) on the 
Galaxy, early devel opment of the method was undertaken by 
Ivan der Kruitl dl98ll) who applied it to low resolution observa- 
tions of NGC891, concluding that the halo was not as flattened 
as the stellar disk. It was then applied to several galaxies in the 



10 



J.C. O'Brien et al.: The dark matter halo shape of edge-on disk galaxies IV 



1990's, most notably the careful study of the very nearby Sc 
galaxy NGC4244, which found a highl y flatte ned halo with by 
q = 0.2^ out to radii of * 2R 25 (Oiling! [19961) . All applications 
of the flaring method h ave indicated highly flattened halo distri- 
butions with q < 0.5 dBecquaertl 1997b iBecauaert & Combesl 
119971; ISickind. 119971) . Recently, iBaneriee & Jod d2008l) mea- 
sured a flattening of q — 0.4 from flaring of the HI layer in M3 1 . 
This assumed a constant HI velocity dispersion with radius; if 
it is allowed to have a modest decline in the outer disk the flat- 
tening can be made less with q more like 0.5 to 0.6. With the 
exclusion of NGC4244, it may therefore be suspected that the 
assumption of a radially constant gas velocity dispersion has led 
to errors in the derived flattening of the halo. 

Measurements of significant flattening using the flaring 
method initially led to the supposition that perhaps the method 
is systematically biased to flattened halos. Our analysis of 
UGC7321 shows that this is not the case: the gas layer flar- 
ing method is just as sensitive to prolate halos as it is to oblate 
ones. Here, we briefly consider the set of q measurements using 
the flaring met hod that have indicated flat halos. The flattening 
for NGC891 (iBecquaert & Combesl 119971) was estimated from 
VLA observations with a low peak signal-to-noise of 13 ((Rupen, 
1199 lb . The low sensitivity could have led to underestimates of 
the gas density and vertical flaring, thus changing the shape of 
both dK- tot/d z and dK 7S /dz , and thereby q. Except in the case 
of NGC4244 (Qllin5 ll996T) . it is unclear what model was used 
for the radial gas velocity dispersion. An assumption of radially- 
constant gas velocity dispersion could easily skew the derived 
halo shape measurement. 

In some cases it is not clear whether the gas self-gravity 
was included in the mass modelling. Additionally, excluding 
NGC4244, all the previous measurements of the halo flattening 
from the gas layer flaring were performed on large Sb-Sc galax- 
ies with maximum rotation speeds v max between 177 and 295 
km s _1 . As the gas layer flaring is inversely proportional to v max , 
the maximum HI flaring of these galaxies is < 1 kpc, making it 
difficult to resolve unless the galaxy is nearer than 5 Mpc. 

The Galactic q mea surement from the gas flaring by 
lOllin g & Merrifieldi (l2000l) is particularly interesting. They were 
unable to fit the halo with a pseudo-isothermal model, unless the 
Solar radius and rotation velocity are significantly less (R < 7.6 
kpc, ©o < 190 km s _1 ) than the standard values. The uncertainty 
associated with these values translates to a large uncertainty of 
4:4 = 0.8 + 0.3. 

UGC7321 is the least massive galaxy for which the halo flat- 
tening has been measured. The derived /?-band face-on central 
surface brightness is 2.5 times fainter than the B-band measure- 
ment of NGC4244, and the total /?-band luminosity is 4.5 times 
fainter than NGC4244, while its gas layer flares to twice the 
height of NGC4244. The very low stellar mass of UGC7321 
made it an ideal candidate for halo modelling with the gas flaring 
method. 

Although there are now a number of different measurements 
of galactic halo flattening, there is no obvious concentration 
around a particular halo shape or any correlation of halo flatten- 
ing with galaxy morphology. Currently the measured q values 
range from 0. 1 to 1 .4. The low q values for the large Sb galax- 
ies, M31, NGC891 and NGC4013, are puzzling as in these cases 
the stellar density distribution may be more spherical than the 
halo density distribution. It seems unlikely that the galactic ha- 
los could exist in the range of shapes measured, unless the frac- 
tions of the constituent dark matte r types vary sig nificantly from 
galaxy to galaxy. Early work by Dubinski ( 1994) found that in- 
cluding baryon infall in n-body halo simulations led to nearly 



axisymmetric halos. Most n-body sim ulations withou t hydrody- 
namics tend to form prolate halos dSellwoodl [2004); however, 
new work by Dubinski (unpublished) has shown that the inclu- 
sion of hydrodynamical modelling generates halos that are more 
spherical. 

We note here an app lication of HI hydrostatics to our Galaxy 
by iKalberla et al.l (120071) . which illustrates the potential power 
of HI hydrostatics to trace the Galactic potential gradient and 
hence the total dark matter distribution in the Galaxy. Kalberla 
et al. adopted an isothermal velocity dispersion for the Galactic 
HI and found several components of dark matter, including the 
usual extended halo with a mass of about 1 .8 x 10 12 M Q , a thick 
self-gravitating disk with a mass of about 2 to 3 x 10 11 M , and 
an outer dark matter ring with a mass of about 2 to 3 x 10 10 M Q . 
Similar studies in other edge-on galaxies may reveal compara- 
ble substructure in the dark matter, including dark matter rings 
which may be left over from accreted and circularized smaller 
galaxies, drawn down into the disk by dynamical friction. As 
we will argue in the next section, it is important to measure the 
structure, rotation and velocity dispersion of the HI in both R 
and z, to ensure such structures are not artifacts of assumptions 
required to apply the hydrostatics. 

Finally we note that since we submitted the original version 
of this paper, a study of the density distrib ution of dark matter 
halo of UGC 7321 by Baneri eeet al.l (12010) appeared, using the 
rotation curve and flaring of the HI layer derived from the same 
data. In this study the fitting was performed with the halo central 
density, core radius and radial exponential density slope as free 
parameters, but with the halo assumed spherical, the stellar M/L 
of the disk fixed and using values for the HI velocity dispersion 
from Gaussian fits to the position-velocity profiles (typically 7 
to 9 km s~'). These authors also conclude from their work that 
the dark matter halo dominates the dynamics of UGC 7321 at 
all radii, but they do not put constraints on the dark matter halo 
flattening. 



6. Conclusions 

In this study we have shown that it is possible to measure the gas 
flaring and HI velocity dispersion via modelling of the HI dis- 
tribution. Using these methods we found that the small late-type 
disk galaxies in our sample show substantial HI flaring, increas- 
ing linearly with radius in the inner disk and exponentially in 
the outer disk. The HI velocity dispersion has a mean value of 
7 km s _1 , but varies from 4.5 to 12 km s _1 . Our HI modelling 
method is also capable of measuring the vertical variation of the 
HI velocity dispersion given additional HI observations. 

UGC7321, a small low surface brightness Sd galaxy in our 
sample, has the most accurate flaring measurements in our sam- 
ple. We were unable to model the observations using a pseudo- 
isothermal halo. By lowering the asymptotic halo rotation to a 
value corresponding to a true isothermal halo model, we found 
that UGC7321 has a spherical halo density distribution of q = 
1.0 + 0.1. Highly prolate halos (q > 1.2) and highly flattened 
halos (q < 0.6) are strongly excluded if our approximation of a 
true isothermal halo is valid. 

Our mass modelling analysis assumed that the HI gas veloc- 
ity dispersion was vertically isothermal, as no measurements of 
the vertical variation of the HI gas velocity dispersion are as yet 
available^ If the HI velocity dispersion is in fact vertically de- 



2 Our measurements of the HI gas velocity dispersion used the ver- 
tically averaged HI distribution, ie. they are the luminosity-weighted 
mean velocity dispersion as a function of radius. 
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clining, this would lead to a larger estimated vertical gradient 
of the total vertical force, which may allow a pseudo-isothermal 
model for the halo. 

UGC7321 is a gas-rich galaxy (MhiILr = 2.2), with a very 
low stellar mass galaxy (M = 3 x 10 s M Q ), four times less mas- 
sive than the gas disk. The /?-band stellar mass-to-light ratio of 
UGC7321 is very low at M/L R < 0.2. Mass modelling of the ver- 
tical force distribution showed that vertical force fitting provides 
a much stronger constraint on the stellar mass-to-light ratio than 
the standard method of radial force fitting via rotation curve de- 
composition. 

Two important assumptions in this work need to be tested 
further. The first is that the HI velocity dispersion is isothermal 
in z. For a definitive estimate of the R and z-components of the 
total potential gradients from HI hydrostatics, it is essential to 
have reliable measurements of the HI density, rotation and veloc- 
ity dispersion as a function of both R and z. It should be possible, 
with additional short spacing ATCA observations supplementing 
our data, to measure the HI velocity dispersion as a function of 
z in ESO274-G001 by modelling the HI XV diagram at varying 
heights above the galactic plane. ESO274-G001 is the closest, 
isolated, southern edge-on galaxy at a distance of 3.4 Mpc. In 
the northern hemisphere, UGC7321 is a prime candidate due to 
its high HI mass, despite its larger distance of 10 Mpc. The large 
HI flaring means that the HI could be measured at a height of 400 
pc for radii from 5-11 kpc, and at 700 pc for radii from 9-11 kpc. 
The dwarf Scd galaxy NGC5023 i s also an excellent candidate , 
given its distance of about 8 Mpc (Ivan der Kruit & SearleL 1982). 
For th is galaxy, early flaring measurements by iBottema et al.l 
(11986 ) found that the gas thickness was constant with radius. 
This is a surprising result, because the v max value for this galaxy 
is only about 80 km s _1 , and large flaring might be expected. It 
would be interesting to measure the radial and vertical variation 
of the gas velocity dispersion, gas flaring, and halo shape using 
beter data, as this galaxy has a similar size, HI brightness and 
total mass as UGC7321. 

The other important test is to determine whether a true 
isothermal halo provides a better model than the pseudo- 
isothermal halo for the dark matter in late-type disk galaxies, 
or whether there are better models than either of these. Our 
analysis of UGC7321 has shown that the vertical gradient of 
the vertical force provides a significantly stronger constraint on 
the halo density distribution than does rotation curve decompo- 
sition. So, this test can in principle be achieved by analysing 
UGC7321 and the other galaxies in our sample with both flat- 
tened pseudo-isothermal and true isothermal halo models. Such 
flattened isothermal halos could be flattened by rotation or by 
anisotropy of the velocity dispersion. This will determine which 
kind of model is better for both the radial halo force as mea- 
sured from the rotation curve, and the vertical force of the halo 
determined from dK-Jdz fitting. 
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